Scattering of low Reynolds number swimmers 



00 

o 
o 

(N 

-t— > 
o 

O 



o 

s 

I 

C 

O 
o 



> 

m 
oo 



oo 
o 



X 



G. P. Alexander, C. M. Pooley, and J. M. Yeomans 
Rudolf Peierls Centre for Theoretical Physics, 1 Keble Road, Oxford, 0X1 3NP, United Kingdom. 

(Dated: October 9, 2008) 

We describe the consequences of time reversal invariance of the Stokes' equations for the hydro- 
dynamic scattering of two low Reynolds number swimmers. For swimmers that are related to each 
other by a time reversal transformation this leads to the striking result that the angle between the 
two swimmers is preserved by the scattering. The result is illustrated for the particular case of a 
linked-sphere model swimmer. For more general pairs of swimmers, not related to each other by 
time reversal, we find hydrodynamic scattering can alter the angle between their trajectories by 
several tens of degrees. For two identical contractile swimmers this can lead to the formation of a 
bound state. 

PACS numbers: 47.63.mf, 47.63.Gd, 47.15.G- 



The motile behaviour of micron sized organisms of- 
fers an insight into a physical environment very differ- 
ent to our own. Micron length scales correspond to low 
Reynolds number conditions where viscous forces domi- 
nate over the effects of inertia. Since Taylor's seminal 
paper [l[ there has been considerable progress in our 
understanding of how low Reynolds number swimmers 
generate their motility @, 0, 0, S @|- In the past few 
years this has included both the development of artificial 
microswimmers Q and a number of simple theoretical 
models d©]. 

A topic of growing interest is the role played by hydro- 
dynamic interactions in determining low Reynolds num- 
ber swimming. These interactions may be expected to be 
substantial because of the long range nature of the fluid 
flow generated by point forces at low Reynolds number, 
and have already been shown to be important in mag- 
netotactic band formation [Io| and in m any aspects of 
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bacterial behaviour near surfaces 
namic interactions between swimmers have been studied 
using a variety of theoretical models, including flagclla 
driven micromachines 
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igella 

rigidly rotating helices , 
squirmers [ 1 dlj . linked sphere swimmers 17 1 and simple 
'body' and 'thruster' models [l§ |. 

A vital concept in understanding the swimming of mi- 
croscopic organisms is that the Stokes' equations, which 
govern zero Reynolds number fluid flows, do not possess 
any intrinsic notion of time. For an incompressible fluid 
of viscosity (i, the fluid velocity u and pressure p satisfy 



^V 2 u - Wp = 







(1) 



and the flow throughout the entire fluid is determined by 
specifying the instantaneous boundary conditions. The 
fluid moves when the boundaries move and stops when 
the boundaries stop. If the motion of the boundaries is 
reversed then the fluid flow is also reversed and each fluid 
element returns to its original position, a phenomenon 
known as kinematic reversibility of Stokes flows. This 
has important consequences for the locomotion of micro- 
scopic organisms, for if their motions are reciprocal, such 



as the opening and closing of a single-hinged scallop 0], 
then kinematic reversibility implies that the forward mo- 
tion of the first half of the stroke is exactly cancelled 
during the second half and there is no net motion, a re- 
sult commonly referred to as the Scallop theorem. 

Kinematic reversibility also implies that when the mo- 
tion of a swimmer (A) is reversed, it produces a second 
swimming stroke (A) just with the swimmer moving in 
the opposite direction. We refer to this as the T-dual 
swimmer, which may be interpreted as the original swim- 
mer going backwards in time. In this Letter we exploit 
the time reversal invariance of the Stokes' equations to 
show that, during any scattering event involving a swim- 
mer and its T-dual, the initial state as t — > — oo is recov- 
ered exactly in the final state as t — > +oo, i.e., the an- 
gle between the swimmers is unchanged by the scattering. 
Experimental verification of the scattering behaviour we 
describe should be possible using biological or fabricated 
microswimmers . 

For simplicity we focus our attention on planar scat- 
tering, however our results naturally generalise to three 
dimensional geometries. Two swimmers, A and B, travel 
along coplanar trajectories with instantaneous swimming 
directions n A (t) and n B (i) (see Fig. [T]). Provided these 
directions are not parallel they generate two straight lines 
which intersect at some point. The angle between these 
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FIG. 1: (Colour online) Schematic diagram of the geometry 
of a planar scattering event: see text for details. 
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lines defines the angle of incidence, a, and the differ- 
ence in the distances, s A , s B , of the two swimmers from 
the intersection point defines the impact parameter, b — 
s A — s B . The limiting values (ct±,b±) = lim t ^± DO (a, b) 
are used to define the initial and final states. The change 
in orientation of an individual swimmer as a result of the 
scattering process is described by the scattering angle 



9 — arccos(n(t — > +00) • n(f — > — 00)) 



(2) 



which we take to be positive if the swimmer rotates in 
an anticlockwise sense and negative otherwise. 

Hydrodynamic scattering may be viewed as providing 
a map from the initial state (a_,6_) to the final state 
(a+,&+). The differences between these two states, de- 
fined by the functions S a :— a+ — a_ and 5^ := b + — b- , 
describe the tendancy for the swimmers to align (8 a < 0) 
or cluster (5b < 0) via hydrodynamic interactions. View- 
ing the entire process backwards in time corresponds 
to the hydrodynamic scattering of the T-dual swimmers 
(B, A) taking the initial state (a+, &+) into the final state 
thereby establishing an isomorphism between 
the scattering of an arbitrary pair of swimmers and the 
scattering of their T-duals. In particular the functions 
§(B,A) are s i m piy related to the functions S^' 3 * 1 



5^ B ' A \a + ,b + ) = -8^ B \a-,b-) 



(3) 



In the case of a pair of mutually T-dual swimmers (B = 
A;A = B) Eq. ([3]) is sufficient to show that (a+,6+) = 
(a_,6_). 

We motivate this result using symmetry arguments. 
During any scattering event the quantity s A + s B changes 
from being large and positive to being large and negative. 
Since it does this continuously it must pass through zero, 
which we use to define the time t = 0. The separation 
between the two swimmers is given by r = (s A +s B )(n A — 
n B )/2 + (s A - s B )(n A + n s )/2, and is orthogonal to the 
direction n A — n B at t = 0. At this instant reversing the 
direction of time, followed by a it rotation about an axis 
parallel to n A — n B and passing through the point mid- 
way between the two swimmers leads to a configuration 
where B, A have the same positions and orientations as 
A, B, respectively. For mutually T-dual swimmers this 
returns the same configuration we started with. It follows 
that A's outgoing trajectory for t > will be given by 
B's ingoing trajectory for t < (with the direction of 
time reversed) and vice-versa, from which we conclude 
that the initial and final states are the same. In addition, 
this construction implies that the swimmers rotate in the 
same direction and with equal scattering angles, 9 A = 
9 B . We therefore refer to these as turn events and show 
an example in Fig. [HI). 

An exception to this scenario occurs if the swimmers 
ever become exactly parallel, n A = n B . However, taking 
t = at this instant and choosing the rotation axis to 
be parallel to n A x r leads to the same conclusion. This 
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FIG. 2: (Colour online) Box: Schematic diagram of the 
Golestanian model swimmer [j| and its T-dual. The swimmer 
is self T-dual if the two arm amplitudes are equal, £ R — £ F . 
Below: Exemplary scattering trajectories of two identical, self 
T-dual Golestanian swimmers obtained using the Oseen ten- 
sor description of the hydrodynamics. In (I) the swimmers 
rotate in the same direction, while in (II) they rotate in op- 
posite directions and exchange trajectories. The rotation axes 
for the symmetry transformations are indicated in the central 
panel, which corresponds to the time t = 0. The initial con- 
ditions were a_ = 30° with &_ = 3.5D for the turn event (I) 
and b- — 2D for the exchange event (II). 



time, since the two swimmers rotate in opposite direc- 
tions, the constraint that a + = q_ can only be met if 
the scattering angles take the values 9 = ±a_ indepen- 
dent of &_. In such an event A will rotate so that its 
outgoing trajectory is parallel to J5's ingoing trajectory 
and vice- versa. We call this an exchange event, an ex- 
ample of which is shown in Fig. [2JII) . Since we expect 
9 —> as &_ — ► co, exchange events can only occur for 
sufficiently small values of 6_. Finally, we comment that, 
since for purely planar scattering there is no way to cross 
smoothly between these two cases, they are necessarily 
separated by some form of discontinuous behaviour. 

There is a subtlety in the foregoing observations. To 
exactly interchange the swimmers as described, the time 
t = must coincide with particular instants during the 
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swimming cycle, otherwise, although the positions and 
orientations of the swimmers will be the same, the stages 
they are at during their swimming strokes will not. We 
have not been able to show generally that the time t = 
does indeed coincide with one of these instances (al- 
though for this not to be the case would imply rather 
peculiar properties for the functions S^ A ' A ^). However, in 
numerical tests using the Golestanian model the swim- 
mers do indeed reach t = at a suitable stage of their 
stroke. Also, since our discussion has mentioned only the 
swimming direction, n, it has been restricted to swim- 
mers that are axisymmetric, requiring only this vector to 
completely specify their orientation. 

These general symmetry considerations provide a 
framework for what can be expected in two body swim- 
mer scattering. In the remainder of this Letter we il- 
lustrate and extend the results by describing the hydro- 
dynamic scattering of a simple model of linked sphere 
swimmers first introduced by Najafi and Golestanian Q. 
Three spheres, each of radius a, are connected by thin 
rods of natural length D, as shown in Fig. [2] By period- 
ically extending and contracting these rods the organism 
is able to swim in the direction of its long axis. Many 
different swimming strokes are possible, but we consider 
here only a simple stroke in which the rods undergo si- 
nusoidal oscillations with amplitudes £, R ,£, F and with a 
phase lag <f> between them. Viewing this motion back- 
wards in time, one sees that T-duality corresponds to 
an interchange of the amplitudes, £, R <-> £ F . Note that 
the swimming stroke is unchanged if the amplitudes are 
equal. This is an example of a special type of swimmer 
which we refer to as self T-dual. Swimmers of this type 
have strokes that 'look the same forwards and backwards 
in time' and yet are not reciprocal: they are time rever- 
sal covariant rather than invariant. A number of simple 
model swimmers commonly referred to in the literature 
are self T-dual, such as Taylor's rotating torus ,1], Pur- 
cell's three link swimmer [3j, the 'pushmepullyou' swim- 
mer Q , a sinusoidally waving sheet [l[ and a rigidly ro- 
tating helical filament 0]. 

The swimming motion of a single Golestanian swimmer 
may be determined analytically using the Oseen tensor 
to describe the hydrodynamics This approach may 
also be applied to determine the interactions between two 
swimmers [13] • These interactions prescribe how the po- 
sitions and orientations of each swimmer are altered by 
the fluid flow generated by the other during one swim- 
ming stroke. By iterating theses changes numerically we 
are able to generate the trajectories of the two swimmers 
during a scattering event. Since this approach is based 
on the Oseen tensor approximation to the hydrodynamics 
it is only valid so long as the separation r between the 
swimmers is large compared to the size of the spheres, 
i.e., a/r <C 1. An important feature of the interactions 
is that they are strongly sensitive to the relative phase 
rj of the two swimmers, which therefore has a significant 
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FIG. 3: (Colour online) Hydrodynamic scattering of two iden- 
tical Golestanian swimmers obtained using the Oseen tensor 
description of the hydrodynamics, (a) The type of scattering 
observed for different values of the initial conditions 
when the swimmers are in phase, rj = 0. (b) Dependence of 
the scattering angle on 6_ for a fixed value of a_ = 30°, cor- 
responding to the dashed line in (a). The change in alignment 
d a is shown as a function of (c) 6_ for a_ = 30°, rj — tt/2 and 
(d) rj for a_ = 30°, b- — 4D, for two extensile (circles) and 
two contractile (crosses) swimmers, (e) Variation of 5 a with 
the difference in amplitudes £ F — £ R for two identical swim- 
mers with q_ = 30°, b- — 4D and rj — tt/2 and (f) exemplary 
trajectory of the bound state formation for £ F — £ R = 0.17ZX 



influence on the type of behaviour that is observed [ljj . 

We consider first two identical swimmers that are in 
phase, rj = 0, and have equal arm amplitudes, £ R — 
£ F , for which the swimming stroke is self T-dual. For 
all trajectories (a+,fe+) = (a_,6_) in accordance with 
the symmetry arguments we have presented above. The 
type of scattering event (exchange or turn) that occurs is 
shown as a function of the two initial conditions a_ , 6_ 
in Fig. OUa), together with a detailed cut showing how 
the scattering angle 6 varies with 6_ for a fixed value of 
a_ = 30° (Fig. EJb)). There is a wide range of initial 
conditions at small values of 6_ for which the scattering 
is of the exchange type and 6 — ±cv_. At larger values 
of &_ the scattering is always of the turn type with the 
scattering angle decaying to zero as 6_ — » oo. In the 
region labelled Collide the swimmers approach so closely 
that the Oseen tensor description of the hydrodynamics 
is no longer valid and we are unable to determine what 
happens during the scattering. In our simulations we 
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took this to occur when the minimum separation between 
any two spheres became less than 10a. 

We now outline how the properties of swimmer scatter- 
ing change when the two swimmers are not mutually T- 
dual and our preceeding symmetry arguments no longer 
apply. In Fig. [3fc) we show the dependence of S a on 
6_ for a_ = 30° and r\ = tt/2. For a pair of iden- 
tical extensile swimmers with (£ , £ F ) = (0.3D, 0.1D) 
the scattering predominantly leads to an increase in the 
angle between the two swimmers, i.e., 8 a > 0. By 
contrast, a pair of identical contractile swimmers with 
(? fl >£ F ) = (0.1D, 0.3D) predominantly shows hydrody- 
namically induced alignment, S a < 0. The relative phase 
is important for these results, as illustrated in Fig. [3jd). 
In particular, the sign of S a is different for < r\ < ir and 
7r < rj < 27r for both extensile and contractile swimmers, 
revealing that there is not a simple, direct link between 
the type of swimming stroke and a tendency for hydro- 
dynamic alignment. Moreover, for 77 = 0, n we find the 
intriguing result that S a — even though the swimmers 
are not mutually T-dual, a result which is also found for 
all other initial conditions (at—, 6_) so long as the swim- 
mers do not collide (not shown). 

One further feature of the scattering of contractile 
swimmers deserves mention, namely the appearance in 
Figs. \3lc) and (d) of a range of scattering events for 
which 5 a takes the constant value —30°, corresponding 
to a + = 0°. This represents the formation of a bound 
state in which the two swimmers are exactly aligned one 
behind the other. An example of swimmer trajectories 
during the formation of this bound state is shown in 
Fig. [2Ie). On either side of the bound state region 8 a 
takes a value somewhat larger than —30°, indicating 
that S a is discontinous at the transition from scattering 
to bound state formation. This qualitative observation 
is supported by the behaviour when the transition is 
approached from a different direction, that of increasing 
the difference in amplitudes £ F — £ R at a fixed value of 
6_ = AD as shown in Fig. [3^d). We vary this difference 
in amplitudes whilst holding the product £ R £ fixed 
to maintain an approximately constant swimming 
speed [H [T3|. As £ F — £ R is increased from zero S a 
decreases smoothly until it reaches the value 0.16Z? 
where there is an abrupt transition to the bound state, 
accompanied by a substantial discontinuity in 6 a . 

We have described the constraints imposed on the hy- 
drodynamic scattering of two swimmers by the time re- 
versal invariance of the Stokes' equations. The most 
striking observation concerns T-dual swimmers, which 
have strokes that map onto each other under time rever- 
sal: for scattering events involving two such swimmers 
the angle and impact parameter between their trajec- 
tories are the same before and after the collision. For 



swimmers unrelated by T-duality we show numerically 
that scattering is complex, with the possibility of changes 
in the angle between the two swimmers of several tens of 
degrees or the formation of bound states. Experiments 
on biological or fabricated microswimmers should show 
these striking differences between pairs of mutually T- 
dual and symmetry unrelated swimmers. In future work 
it may be possible to extend our ideas to more than two 
swimmers, thereby constraining the properties expected 
of large groups or swarms, and providing a connection to 
continuum models. 

We are grateful to Mike Cates, Scott Edwards, Davide 
Marenduzzo, and Vic Putz for enlightening discussions 
about this work. 
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